The X-ray Spectrum of SAX J1808.4-3658 



W. A. Heindl 

Center for Astrophysics and Space Sciences, Code 0424, University of California, San 

Diego, La Jolla, CA 92093 

D. M. Smith 

Space Sciences Laboratory, University of Cahfornia, Berkeley, Berkeley, CA 94720 

ABSTRACT 

We report on the X-ray spectrum of the 401 Hz X-ray pulsar and type I 
burst source SAX J1808. 4-3658 during its 1998 April/May hard outburst. 
The observations were made with RXTE over a period of three weeks. The 
spectrum is well-described by a power law with photon index 1.86 ± 0.01 that 
is exponentially cut off at high energies. Excess soft emission above the power 
law is present as well as a weak Fe-K line. This is the first truly simultaneous 
broad-band (2.5-250 keV) spectrum of a type I burst source in the hard state. 
The spectrum is consistent with other hard state burster spectra which cover 
either only the soft (1-20 keV) or hard ( ^ 20keV) bands, or cover both, but 
not simultaneously. The cut-off power law spectrum resembles that of black 
hole candidates (BHCs) in their low states, observed with RXTE. We compare 
the SAX J1808. 4-3658 spectrum to three BHCs and find that the power law 
is somewhat softer. This suggests that the photon index may provide a way 
to distinguish between low state emission from Galactic black holes and type I 
bursters. 



Subject headings: stars: individual (SAX J1808. 4-3658) — stars: neutron — 
X-rays: stars 



1. Introduction 

SAX J1808. 4-3658 is the first object thought to display both type I X-ray bursts and 
coherent X-ray pulsations. Its low implied magnetic field (B 2 x 10® G, Wijnands & 
van der Klis 1998b) and high spin frequency may make SAX J1808. 4-3658 a missing link 
in the evolution of the millisecond radio pulsars. It was discovered in observations of the 
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Galactic center region made during 1996 September 12-17 with the Wide Field Camera 
(WFC) on BeppoSAX ( |in't Zand et al. 1998| ). During six days of observations, the flux 
level was ~ 50-100 mCrab (2-10 keV). Earlier and later observations which did not detect 
the source limit the outburst duration to between 6 and 40 days. This duration was 
conflrmed with data from the All Sky Monitor (ASM) on the Rossi X-ray Timing Explorer 
{RXTE}, which detected SAX J1808.4-3658 for about 20 days beginning 1996 September 8 
( |in't Zand et al. 1998|) . Two type I X-ray bursts were also detected during the BeppoSAX 



observations, making the identiflcation of the source as a neutron star in a low mass X-ray 
binary highly probable. Assuming that these bursts reached the Eddington luminosity for a 
1.4M0neutron star implies a distance of ~4kpc. 

Following the 1996 outburst, SAX J1808. 4-3658 remained undetected until a slew of the 
RXTE pointed instruments on 9 April 1998 serendipitously detected a source (designated 
XTE J1808-3658) whose location is consistent with the BeppoSAX error region ( IMarshalll 



1998|) . The flux level at this time was ~50mCrab (2-10 keV), corresponding to a luminosity 
of 1.5 X lO^^ergs/s at a distance of 4kpc. Twenty-one RXTE pointed observations over the 
next 4 weeks saw the flux increase to 60mCrab (2.5-20 keV) and decrease approximately 
exponentially with a time constant of about 10 days (see Figure jT]). After 26 May, the 
source dimmed rapidly by a factor of ~5 in 2 days . 

Timing analyses of i^XTE'/Proportional Counter Array (PCA) data from 11 April 
1998 revealed that SAX J1808. 4-3658 is an X-ray pulsar with a frequency of 401 Hz, 
making it the flrst accretion-powered millisecond X-ray pulsar ([Wijnands fc van der Klis 



19984 [Wijnands fc van der Klis 1998b| ). The pulsed amplitude was quite low, only ~4% 







RMS (2-60 keV). Chakrabarty & Morgan (1998a), also using PCA data, detected the 
binary orbit. They derived an orbital period of 7249.119(1) s, a projected semimajor axis 
of axSin^ = 62.809(1) It-ms, and an X-ray mass function of 3.85 x 10~^ Mq ( [Chakrabarty 



fc Morgan 1998 b| ). They also placed upper limits on the pulse frequency derivative and 
the eccentricity (< 5 x 10^^) of the orbit. The very small mass function implies that the 



companion mass is < 0.18 Mq for a neutron star less massive than 2 Mq ( Chakrabarty & 
Morgan 1998b[ ). 



During the recent X-ray outburst, optical imaging of the SAX J1808. 4-3658 fleld 
revealed an object with magnitude V = 16.6 that was not present in the Digitized Sky 
Survey (DSS) to a limiting magnitude of V > 19 ( [Roche et al. 199^ ). This object was 
conflrmed as the likely optical counterpart of SAX J1808. 4-3658 when multiple V-band 
exposures covering the 2hr binary orbit showed "roughly sinusoidal" variability of 0.12 mag 
( [Giles, Hill fc Greenhill 19981) . 



Early work with the RXTE data indicated that the SAX J1808. 4-3658 spectrum was 
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Fig. 1.— The light curve of the 1998 April outburst of SAX J1808.4-3658 in 
three energy bands from the ASM, PCA, and HEXTE on RXTE. The ASM 
data are the publicly available daily averages provided by the ASM/RXTE team 
( |http:/ /heasarc.gsfc.nasa.gov/docs/xte/asm_products.html| ). 
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Crab-like and continued unbroken to energies greater than 100 keV ( Jieindi, Marsden fc 



Blanco 1998| , piles, Hill fc Greeniiill 1998| . In this letter, we perform detailed spectral 
studies and show that the spectrum is somewhat harder than the Crab at low energies 



< 



30 keV) and is exponentially cut off at higher energies. 



2. Observations and Analysis 

The recent outburst was the subject of a monitoring campaign with the PCA and 
the High Energy X-ray Timing Experiment (HEXTE) on RXTE. The PCA ( |Jahoda et al 



1996D is a set of 5 Xenon proportional counters (2-60 keV) with a total area of ~7000cm^. 



The HEXTE consists of two clusters of 4 NaI(Tl)/CsI(Na) phoswich scintillation counters 
(15-250 keV) totaling -1600 cm^ area QGruber et al. 1996| , [Rothschild et al. 1998|) . The two 



clusters alternate rocking between source and background fields to measure the background. 
The PCA and HEXTE share a common 1° full width half maximum field of view. In this 
Letter we discuss the X-ray spectra obtained during 13 observations made between 1998 
April 11 and 25. Figure ^ shows the flux history of SAX J1808. 4-3658 in three energy 
bands. The outburst peaked between April 10 and 12 and declined over the next 20 days. 

For each observation, we accumulated PCA and HEXTE spectra without regard to the 
binary orbit and the X-ray pulse phase. These spectra are quite hard, being superficially 
similar to the Crab Nebula and pulsar (a power law with photon index of ~2) and 
extending to over 100 keV ([Heindl, Marsden fc Blanco 1998| , [Gilfanov et al. 1998| ). The 



PCA background was estimated using PCABACKEST version 1.5, with the background 
model based on blank-sky pointings. We then used XSPEC to fit various models to the 
observed counts spectra. PCABACKEST and XSPEC are standard NASA software tools. 
In all fits the relative normalizations of the PCA and the two HEXTE clusters were taken 
as free parameters, owing to uncertainties ( ^ 5%) in the HEXTE deadtime measurement. 
We then verified that the fitted relative normalizations were consistent with those derived 
from fits to the Crab. 

No single component model could adequately fit the full 2.5-250 keV spectrum. In 
particular, excess emission above a power law at low energies ( ^ lOkeV), Fe-K line 
emission, and (for the more significant observations) a cutoff at high energies ( ^ 35keV), 
are all required. Several complex models fit the spectra with acceptable x^. For example, 
the broad-band continuum (above the Fe-K line) could be fit by any of a broken power 
law, a power law with an exponential cutoff above ~30keV, or a Comptonized spectrum 
( ^unyaev fc Titarchuk 19801) . Further, a black body, a disk (multicolor) black body (see e.g. 
Mitsuda et al. 1984 ), and a thermal bremsstrahlung model were all adequate to reproduce 
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the low energy excess. 

Because it provided the best fits and is simple in form, we concentrated on a model 
made up of a power law with an exponential cutoff at high energies {E~^, E < Ecuu 
E~^ X e"(^"-^™*)/'^-f°", E > Ecut) plus a Gaussian Fe-K hue and a disk black body. Low 
energy absorption due to intervening interstellar and/or local gas was also allowed. 

Fits to the individual pointings showed no evidence for spectral variability prior to the 
rapid dimming of the source which began around MJD 50929. Spectra after this time may 
have been slightly softer, but the change was not highly significant. We therefore added all 
the data prior to that date to obtain the most significant spectrum. The PGA data below 
25keV are of extremely high statistical significance, so that uncertainties are dominated 
by small systematics in the instrument response matrix. We applied systematic errors to 
the spectrum of between 0.5 and 2% of the count rate per channel, inferred from fits to the 
Grab. In addition, we fixed the Gaussian width of the Fe-K line at 0.1 keV (a), narrow 
compared to the PGA energy resolution. 

3. Results 

The total spectrum for observations between MJD 50914 and 50928 together with the 
best fit model and the inferred incident spectrum of the form described above are shown in 
Figure |^. The model parameters are listed in Table |I]. Although we included systematic 
uncertainties in the PGA data, we were concerned that the residuals to a simple power 
law model may have been caused by small errors in the response matrix. So, to verify 
our results, we divided the PGA data by a Grab Nebula and Pulsar spectrum (provided 
by K. Jahoda of the PGA team) which was obtained on MJD 50914. The ratio is shown 
in Figure |. Because the Grab is a very bright source, the statistics of its spectrum are 
negligible in the ratio. Since the spectrum of SAX J1808. 4-3658 is similar in overall shape, 
the ratio gives a fairly matrix-independent picture of features which differ from the Grab. 
In particular, the overall positive slope is due to the harder power law in SAX J1808. 4-3658, 
and the soft excess and an iron line are clear below lOkeV. While we remain cautious of the 
physical interpretation of the soft excess as a disk black body and the quantitative iron line 
model parameters, we are confident that both the soft excess and the iron line are present. 

Figure |^ also shows that the resemblance to the Grab is at best superficial. In contrast 
to the power law fits to the individual observations made by Gilfanov et al. (1998), we find 
strong deviations from such a simple model (Xred = 22.4 with 94 degrees of freedom (DOF)), 
even in the limited range of 3-100 keV which they used. The most important excursions are 
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Table 1: Spectral fits to the X-ray burster SAX J1808. 4-3658 and three black hole candidates. 





SAXJ1808.4-3658 


Cyg X-1 


IE 1740.7-2942 


GRS 1758-258 


Nh (xl022cm-2) 


n 07+0.26 
U-J ' -0.21 


0.3 ±0.2 


8.6 ±0.6 


1.3 ±0.3 


TinQikeV) 


n qq+0.08 

'j-yy-0.02 


1.4 ±0.1 


1.15 ±0.30 


1.4 ±0.2 


KdbbB 


13+2 


18 ±4 


1 O+0.8 
-■-■^-0.5 


0.8711^ 


Epe (keV) 




-0.09 


5.7 ±0.70 


6.281°:^?'= 


EWpeQ (eV) 








77 ±9"= 


rg 


1.86 ±0.01 


1.488 ±0.014 


1.53 ±0.06 


1.54 ±0.04 


Ecut (keV) 


34 ±4 


34 ±5 


19.1 ±3.5 




Efoid (keV) 




230 ± 20 


116 ±20 


185;^^ 


Ll~20kevQ 


2.5 


5.2 


7.0 


5.7 


L20-200kev'' 


1.9 


11 


15 


12 



"Temperature at the disk inner edge. 

''Normalization of the disk blackbody component: Kdbb — {Rin/ Diq)cos9. i?i„ is the disk inner radius in 
km, Dio is the distance in units of lOkpc, and 6 is the angle of the disk. 

'^Additional systematic uncertainties apply due to the subtraction of diffuse Galactic plane Fe-K emission. 

'^Fe-K line equivalent width. 

"^Power law photon index. 

■''Luminosity (10^'' ergss^^). Assuming distances of SAX J1808. 4-3658: 4kpc, Cyg X-1: 1.9 kpc, 
IE 1740.7-2942: 8.5 kpc, and GRS 1758-258: 8.5 kpc. 
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Fig. 2.— The spectrum of SAX J1808.4-3658 between MJD 50914 and 50928. The upper 
panel shows the measured and best fit counts spectra (data points and histograms) and the 
inferred incident photon spectrum (smooth curve). 
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at low and high energies where we find a soft excess and an exponential cutoff respectively. 
The fit in Table |1] for SAX J1808.4-3658 had = 1-65 for 99 DOF (2.5-250 keV). If we 
disallow the exponential cutoff and refit all the other parameters, we find Xi-cd ~ 4-34 for 
101 DOF. Applying an F-test ( |Bevington 1969 ) to the ratio of the improvement in and 



of the full model (divided by 2 and 99 DOF, respectively), we find that the probability 
that the cutoff is unnecessary is only 7 x 10~^^. 

We also analyzed the PGA spectra in eight bins by orbital phase. Due to the limited 
amount of data, each phase bin samples the different parts of the decaying lightcurve to 
different degrees. To compensate for this effect, we fit an exponential decay to the lightcurve 
and looked for additional variations. We find an apparent intensity variation with orbit, 
in agreement with Chakrabarty & Morgan (1998b): roughly sinusoidal, with an amplitude 
of 2%, and a minimum when the neutron star is behind its companion. The modulation 
is several times larger than could be explained by uncertainties in the background model. 
We find no differences in the spectrum with phase: the disk blackbody and power law 
components share the extinction equally, and the absorption column is unchanged. This 
may imply that the modulation is due to Compton scattering in a thin, ionized intervening 
medium (e.g. the ablated wind suggested by Chakrabarty & Morgan (1998b)), which 
would scatter photons out of the line of sight nearly independently of energy. Alternately, 
the modulation could be an artifact of the limited amount of data, caused by random 
correlations between the orbital phase and small variations in the accretion rate: the total 
amount of data available is only about 14 orbits. 



4. Discussion 

4.1. Comparison to Type I Burst Sources 

Prior to this work, the only bursters with nearly simultaneous coverage of the X-ray 
and hard X-ray bands during episodes of hard outburst were Cen X-4 and 4U 1608-52. 
In 1979, Cen X-4 was observed with Hakucho (1.5-30 keV) ( [Matsuoka et al. 1980| ), Ariel 
5 (3-6 keV) ( |Kaluzienski,Holt fc Swank 19801) , and Prognoz 7(13-163 keV) ( [Bouchacourt 



|et al. 198^ ). Bouchacourt et al. (1984) fit the Prognoz 7 spectrum to a hard power law 
(photon index fixed: P = 1) times an exponential cutoff with a folding energy of ~50 keV 
- not unlike the RXTE SAX J1808. 4-3658 spectrum. They also noted that the low energy 
extrapolation of the Prognoz 7 data fell below the simultaneous Ariel 5 spectrum, suggesting 
soft excess emission, although this could be due to the true power law differing from E~^. 
Zhang et al. (1996a) used the Ginga spectrum ([Yoshida et al. 1993| ) from the middle 



of a ~170 day (1991 June - December) outburst to constrain the low energy portion of 
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the 4U 1608-52 spectrum observed by BATSE. The Ginga data showed a power law with 
photon index of —1.75 ± 0.01, and together with the BATSE data the spectrum was fit by 
a Comptonized model ( |Sunyaev fc 'I'itarchuk 195I1| ) with kT = 23keV and optical depth 
4.4 (for a spherical geometry). For SAX J1808. 4-3658, we find nearly identical values of 
22.0lo;gkeV and 4.02^0^5 for the electron temperature and optical depth. Thus, the broad 
band spectra of all three well-studied hard state burst sources are quite similar. 

Hard X-ray emission from other bursters has been observed (without concurrent soft 
X-ray coverage) with BATSE and SIGMA (see Barret, McClintock, & Grindlay 1996 and 
references therein). When these data are fit to broken power law spectra, photon indices of 
between 2.5 and 3 above the break are generally found. Fitting the HEXTE data above 
40keV gives an index of 2.4 ± 0.1 for SAX J1808. 4-3658, in reasonable agreement with the 
high energy spectra of other low state bursters. 

4.2. Comparison to Black Hole Candidates 

Historically, hard emission extending to ^ 100 keV from X-ray binaries has been 
considered a "black hole signature", distinguishing black hole sources from low magnetic 
field neutron stars. Clearly with several burst sources emitting hard X-rays, the situation 
is not so simple. However, it may still be possible to separate these neutron stars by other 
properties of their emission. For example. Barret, McClintock, & Grindlay (1996) noted 
that neutron stars only produce hard tails when the tails dominate the emission, as opposed 
to black holes, which can produce hard tails along with very bright ultrasoft emission. They 
also noted that the luminosity when the hard emission dominates can be higher for black 
holes than for neutron stars. SAX J1808. 4-3658 is consistent with a neutron star by both 
these criteria, showing only a very weak soft component and having a modest luminosity 
(~ 5 X 10^^ ergs s~^ assuming a distance of 4kpc). Given that distances are often poorly 
known and that black hole candidates can also be seen at low luminosities, luminosity alone 
is an insufficient means to distinguish between source types. 

To see if the spectrum alone can serve the purpose, we fit RXTE observations 
of the persistent low state black hole candidates (BHCs) Cyg X-1 (1998 April), 
IE 1740.7-2942(1996 March), and GRS 1758-258 (1996 August) to the same model used 
for SAX J1808. 4-3658. These fits are listed in Table |I]. One feature of the BHC spectra 
clearly differs from the bursters: the photon indices are significantly harder in the BHCs. 
The BHCs power laws are between 1.4 - 1.6, while the SAX J1808. 4-3658 index is 1.86. 
From Ginga observations, 4U 1608-522 had an index of 1.75 during the period of hard 
emission described above. Thus, from this limited sample, the broad-band spectrum 
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separates the BHCs and neutron stars via the hardness of the power law. Given the 
reahties of intcrcalibrating multiple instruments and the modest separation of the indices, 
this characteristic may be difficult to use for comparisons between different observatories. 
However, it will clearly be useful for observations made with a single instrument complement 
or between instruments which are well intercalibrated. 



5. Conclusions 

The spectrum of the unique 401 Hz pulsar and type I burster SAX J1808. 4-3658 
during its recent outburst was quite hard. The photon spectral index was 1.86, with a slow 
cutoff at high energies. There was clear evidence of excess soft emission and a weak Fe-K 
line. These RXTE observations have provided the highest quality broad-band spectrum 
of a suspected X-ray burster during a period of hard emission. The spectrum is in good 
agreement with previous observations of other bursters made over more limited energy 
bands and with lower statistical significance. By comparing to observations of low state 
BHCs, we find that BHCs and bursters can be distinguished by the slope of their power law 
emission. In particular, the photon indices of bursters are greater (i.e. softer) by about 0.3, 
even though the overall spectral shapes are quite similar. 

Thanks to Keith Jahoda for his help with the PCA response matrix and for providing 
the Crab data. We would also like to acknowledge the RXTE Science Operations Center for 
scheduling these Target of Opportunity observations. This work was supported by NASA 
grant NAS5-30720. 
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